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Abstract. We discuss the possibility that the stellar sys- 
tem Cyg OB2 No. 5 can be a gamma-ray source in the light 
of recent EGRET and radio data. This system is formed 
by an 07 la + Ofpe/WN9 contact binary. A third star, 
probably a BO V star also associated with the system, is 
located at ~ 1700 AU from the primary. We estimate the 
expected gamma-ray luminosity from the colliding winds 
region, the terminal shock of the wind, and the unstable 
zone at the base of the wind, and conclude that, under 
very reasonable assumptions, Cyg OB2 No. 5 can gener- 
ate about a half of the gamma-ray flux detected from the 
positionally coincident source 3EG J2033+4118. We sug- 
gest, then, that other O stars belonging to the association, 
also placed within the 95 % probability EGRET location 
contour, could contribute to the observed gamma-ray flux. 

Key words: Gamma rays: theory - Gamma rays: obser- 
vations - Stars: early-type - Stars: winds, outflows 



1. Introduction 

High-energy emission from early-type stars has been often 
discussed from a theoretical point of view (e.g. Casse & 
Paul 1980, Volk & Forman 1982, White 1985, Chen & 
White 1991, White & Chen 1992, Eichler & Usov 1993). 
The observational evidence for its existence, however, has 
remained rather elusive till now. 

The recent publication of the Third EGRET (3EG) 
Catalog by Hartman et al. (1999) provides new infor- 
mation that can be very useful for the identification of 
possible gamma-ray emitting stellar systems. There are 
170 unidentified sources in the new catalog, and approx- 
imately half of them are concentrated near the galactic 
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plane (e.g. Gchrcls ct al. 2000). Romero, Benaglia & Tor- 
res (1999) have studied the positional coincidence between 
3EG sources and different types of galactic objects such as 
supernova remnants (SNRs), OB star associations (usu- 
ally considered as tracers of pulsar concentrations), and 
massive stars like Wolf-Rayet (WR) and Of stars, which 
have very strong stellar winds. They have found that 10 
gamma-ray sources in the 3EG catalog have WR or Of 
stars within their 95 % confidence location contours. The 
probability of this being a mere effect of chance is in the 
range 10 -2 — 10~ 3 . Although such probabilities are not as 
low as those found for SNRs and OB associations, they are 
suggestive enough as to encourage an a posteriori analy- 
sis of some interesting candidates. In particular, Romero 
et al. (1999) mentioned the possibility that the system 
Cyg OB2 No. 5 could be physically linked to the source 
3EG J2033+4118, and hence that it could be the first stel- 
lar system to be detected at gamma-rays, involving other 
than WR stars. In this paper we discuss whether Cyg OB2 
No. 5 can actually generate the observed gamma-ray flux 
and which are the relevant radiation mechanisms for this 
system at high energies. 

Our results will be relevant not only for this particu- 
lar system but also for other interacting binaries like WR 
140 (Romero et al. 1999, Zhekov & Skinner 2000), which 
is also positionally coincident with a gamma-ray source, 
and WR 147, for which there are detailed radio observa- 
tions (Contreras & Rodriguez 1999). The interest of stellar 
systems as potential gamma-ray sources has been recently 
reinforced by the suggested association of 3EG J1824-1514 
with an X-ray binary by Paredes et al. (2000). 

We shall begin our study by reviewing the main char- 
acteristics of Cyg OB2 No. 5. 

2. The stellar system Cyg OB2 No. 5 

Cyg OB2 No. 5 (also known as V729 Cyg, BD 
+40°4220) is an evolved contact binary system with a pe- 
riod of 6.6 days (Hall 1974) located at a distance of ~ 1.8 
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kpc (e.g. Waldron ct al. 1998). The primary star in the 
system is a very massive star classified by Walborn (1973) 
as 07 Ianfp, although it should be taken into account that 
a sizable fraction of the line emission arises in the inter- 
action region between the two stars and is not intrinsic 
to the stars themselves. Most of the spectral features that 
justify the "nfp" tag are actually emission lines that form, 
to a large extent, in the wind interaction region of the 
close binary. 

Bohannan and Conti (1976) also classified the sec- 
ondary star as an Of star with similar luminosity, since 
each star seems to contribute to the He II A4686 emis- 
sion. This result has been recently reconsidered by Rauw 
et al. (1999), who have presented a detailed analysis of an 
extensive set of spectroscopic observations finding a new 
orbital solution. These authors showed that the observed 
properties can be explained if the secondary is a transition 
Ofpe/WN9 star, with a much slower terminal wind veloc- 
ity than the primary. Consequently, at some distance from 
the system, the wind properties are basically determined 
by the primary, with a terminal velocity of Uqo ~ 2200 km 
s" 1 (Bieging et al. 1989). 

Cyg OB2 No. 5 is well known as a rather powerful 
and variable radio source (Abbott et al. 1981; Persi et al. 
1985, 1990; Miralles et al. 1994). The radio source seems 
to switch between states of high and low emission with a 
period of about 7 years. Miralles et al. (1994) detected, 
in addition to a strong source coincident with the contact 
binary, a weak radio component, possibly nonthermal, lo- 
cated at ~ 1 arcsecond towards the NE of the binary sys- 
tem. They proposed that the radio variability is produced 
by relativistic particles accelerated in the wind. The weak 
radio component, however, is too far from the binary sys- 
tem as to be the radio counterpart of the secondary. In its 
approximate position, a star of magnitude 13-14 has been 
previously reported by Herbig (1967). 

Recently, Contreras et al. (1997) have obtained accu- 
rate optical (HIPPARCOS and Roque de los Muchachos) 
and radio (VLA) measurements that clearly show that the 
weak radio source is nonthermal and does not coincide 
with the optical position of the star detected by Herbig 
(1967). This star, which seems to be a BO V star, is a 
third component of the system, with a long period. The 
fact that the weak radio source is located between the 
contact binary and the BO V star, closer to the latter, 
was interpreted by Contreras et al. (1997) as the effect of 
the particle acceleration in the region of colliding winds. 
In this region, electrons locally accelerated radiate syn- 
chrotron emission, which is detected at centimeter wave- 
lengths as the weak radio component. 

The system Cyg OB2 No. 5, then, seems to be formed 
by three stars: an 07 la + Ofpe/WN9 contact binary and 
a BO V star in a larger orbit. When evaluating the gamma- 
ray production in Cyg OB2 No. 5, we shall consider the 
contact binary, whose winds are dominated by the 07 la 



star, as the primary object in the presumed gamma-ray 
emitting system, and the BO V star as the secondary. 

In Table 1 we summarize the stellar parameters 
adopted in this paper. The mass loss rate of the primary 
has been derived from Ha, infrared and radio observa- 
tions, giving values in the range 1.9 to 2.6 x 10~ 5 M©yr _1 
(e.g. Leitherer ct al. 1982, Persi et al. 1990). Contreras et 
al. (1997) propose an extreme M = 5.5 x 10~ 5 M Q yr _1 
that should be considered rather as an upper limit. In 
our present work we adopt an average value of 2.5 x 10~ 5 
MQyr -1 . The mass loss rate for the secondary results from 
that adopted for the primary and geometrical consider- 
ations on the location of the colliding wind region (see 
Section 4.1). The effective temperatures and the stellar 
luminosities are taken from Vacca et al. (1996). There 
are not direct measurements of the stellar magnetic fields, 
although they arc usually thought to be of the order of 
a few Gauss (e.g. Maheswaran & Cassinclli 1992). We 
have adopted values of 2 and 1 G for the primary and 
the secondary, respectively. A larger field would produce 
a turnover in the synchrotron emission of the primary at 
a few GHz that is not observed, whereas a smaller value 
would shift this turnover to hundreds of GHz, resulting in 
an excess of non-thermal emission at high frequencies. Our 
assumptions are in agreement with White & Chen (1992) 
and allow a fair description of the observed radio spec- 
trum for the primary in terms of White's (1985) model 
(see Section 4.3). 

3. The gamma-ray source 3EG J2033+4118 

According to the 3EG catalog (Hartman ct al. 1999), the 
best estimated position of 3EG J2033+4118 is at (I, b) w 
(80.27°, 0.73°), i.e at 0.27 degrees from Cyg OB2 No. 5 
(see Figure 1). The summed EGRET flux (E > 100 McV) 
for cycles from 1 to 4 (i.e. for the time span 1991-1995) 
is (73.0 ± 6.7) x 10~ 8 ph cm~ 2 s" 1 . If we assume that 
the source is at the distance of Cyg OB2 No. 5 (1.8 kpc), 
an integrated gamma-ray luminosity of ~ 2.4 x 10 35 erg 
s" 1 is implied for the EGRET energy range (100 MeV - 
20 GeV). The photon spectral index is T = 1.96 ± 0.10 
(N(E) oc E~ r ). For all calculations we assume r = 2. 

In order to estimate the degree of variability of the 
gamma-ray emission we have adopted the method used 
by Zhang et al. (2000), Romero et al. (2000) and Torres 
et al. (2000). A variability index / is defined such that 
/ = /i s / < n > p , where /i s = 100 a < F > _1 is the fluc- 
tuation index of the gamma-ray source and < fi > p is 
the averaged fluctuation index for all known gamma-ray 
pulsars, which arc usually considered as a non- variable 
population. F is the flux and a its standard deviation. A 
variable source, then, should satisfy that J > 1. For 3EG 
J2033+4118 we obtain / = 1.3, so, taken into account the 
errors involved (la = 0.5), the source cannot be classified 
as clearly variable, at least over the timescale of the ob- 
servations (~ 4 yr). This do not preclude that the source 
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could experience variability on longer timescales, for in- 
stance with a period similar to what is observed at radio 
frequencies in Cyg OB No. 5. (i.e., <~ 7 yr). We shall have 
to wait until the next generation of gamma-ray satellites 
to get a confident evaluation of the temporal behaviour of 
this source. 



Fig. 1. Location map of the gamma-ray source 3EG 
J2033+4118. The contours represent 50%, 68%, 95% and 98% 
statistical probability that the source lies within that contour. 
The position of the stellar system Cyg OB2 No. 5 is indicated 
with a star symbol 



4. Gamma-ray production in Cyg OB2 No. 5 

The main mechanisms that can generate gamma rays at 
energies E > 100 McV in a stellar environment are inverse 
Compton (IC) scattering, relativistic bremsstrahlung, and 
7T°-decay from hadronic interactions. All these mecha- 
nisms require the existence of a relativistic population of 
charged particles (electrons or positrons in the first two 
cases, protons or ions in the latter). Effective particle ac- 
celeration up to high energies requires, in turn, the medi- 
ation of strong shock fronts (e.g. Bell 1978a,b; Blandford 
and Ostriker 1978; Blandford 1980; Drury 1983): the parti- 
cles are scattered by magnetic irregularities carried by the 
converging flow on either side of the shock wave, leading 
to a first-order Fermi process that results in a power-law 
spectrum of relativistic particles. 

In a system such as Cyg OB2 No. 5 there are several 
places where strong shocks can be formed. In particular, 
strong shocks are expected in the colliding winds region 



located at some point between the primary and secondary 
stars (i.e. between the contact binary and the BO V star), 
at the point where the momentum flux of the primary (and 
stronger) wind decline to the value of interstellar pressure 
(i.e. the terminal shock), and at the base of the wind where 
line-driven instabilities are thought to occur. These shocks 
locally re-accelerate the particles making them energetic 
enough as to generate the gamma-rays. We shall discuss 
below in major detail the possible gamma-ray production 
in each of the mentioned regions. 

4-1- Region of colliding winds 

The supersonic winds from the primary and secondary 
stars flow nearly radially out to the shocks. The geome- 
try of the region where the winds collide is described by 
Eichler & Usov (1993) in terms of the star separation D 
and the parameter -q = Miv^^l 'M\Vao,\ as: 



r\ = 



1 + r? l/2 



D. 



T2 = 



n 



1/2 



1+^/2 



D. 
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Here, a subscript "1" stands for the primary whereas the 
subscript "2" denotes the secondary. M is the mass loss 
rate, v oo is the wind terminal velocity, and r the distance 
from the star to the colliding winds region. Since rj < 1, 
this region will be closer to the secondary. 

Eichler & Usov (1993) have shown that particles lo- 
cally accelerated at the strong shock front formed in the 
wind collision can generate observable synchrotron radio 
emission. Contreras et al. (1997) have identified the weak 
non-thermal radio source in Cyg OB2 No. 5 with this 
emission. The same population of relativistic electrons 
that emits the synchrotron radiation should produce high- 
energy gamma-rays through IC interactions with the UV 
photons from the secondary. The frequency of the IC pho- 
tons is given by vie = 4/37 2 ^*, where 7 is the Lorentz 
factor of the relativistic electrons (7 ~ 3.5 x 10 3 for elec- 
trons that produce synchrotron photons of 5 GHz in a field 
of <~ 10~ 4 G) and is the frequency of the seed stellar 
photons. The ratio of synchrotron to IC luminosity is in 
this case (e.g. White & Chen 1995): 



^IC 



840- 



(2) 



where B s is the magnetic field in Gauss, r 2 the distance 
to the secondary in AU, and L 2 is its luminosity in Lq 
units. 

From geometric considerations, Contreras et al. (1997) 
derive rj ~ 0.04, and then r 2 ~ 300 AU. The inte- 
grated synchrotron luminosity between the two observ- 
ing frequencies of the nonthermal radio component is 
L syn ~ 4.5 x 10 27 erg s _1 . Taken into account that for 
a B0 V star the peak of the photon distribution is at 

~ 2 x 10 15 Hz, we have that only electrons with Lorentz 
factors in the range 3.0 x 10 3 - 4.3 x 10 4 will contribute 
to the IC gamma-ray flux in the EGRET energy range. 
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In order to estimate the total synchrotron luminosity of 
these particles, we must first evaluate the magnetic field 
at the colliding winds region. We shall assume that the 
external magnetic field of the star in the absence of stellar 
wind is dipole, and that in the presence of wind it obeys 
the standard B oc radial dependence for large r given 
by Eichler & Usov (1993): 



B » £?*- 



Kot R 



valid for r > i?* — — . 

TA T Vrot 



(3) 



where J?* is the stellar radius, £?* the surface magnetic 
field, V TO t the surface rotation velocity, and ta ~ R* 
the Alfven radius. Using the parameters listed in Table 
1 we find that in the colliding winds region the field is 
B <~ 2 x 10 -5 G, a value about a factor 10 larger than 
the average magnetic field in the standard ISM (Spitzer 
1998). We shall adopt, taking into account the compres- 
sion produced by the strong shock formed in the wind 
collision, a field B s = 10~ 4 G. Then, the integrated syn- 
chrotron luminosity of the electrons that contribute to the 
IC gamma-ray luminosity is L syn « 5 x 10 29 erg s _1 . From 
Eq. (2), we obtain that the corresponding inverse Comp- 
ton luminosity with seed photons from the secondary in 
the EGRET energy range is ~ 5 x 10 34 erg s _1 . If we also 
consider the contribution from the photon field of the pri- 
mary, we get a total IC luminosity of the colliding wind 
region of ~ 8 x 10 34 erg s _1 . 

Additional contributions from bremsstrahlung and 
hadron interactions can be expected from the colliding 
winds region. The density of particles in this region is 
(Eichler & Usov 1993): 



Mi 



47rrfw 00i iTO p /x 

3 x 10 9 / Mi 
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Here, m p is the proton mass and jj, is the mean molecular 
weight of the gas, taken as 1.3. Using Eq. (4) and assuming 
again a compression factor of 4 for the post-shock region, 
we estimate a number density n ~ 2.9 x 10 3 cm~ 3 . This 
means that the total mass in the colliding winds region 
is ~ 2.2 x 1CU 6 Mq. In this calculation we have assumed 
a spherical geometry with a radius equal to the average 
value of the semi-axes in Contreras et al.'s (1997) 6-cm 
elliptical image of the secondary radio source. 

The contribution from 7r°-decays, consequently, should 
be very low, even if an important enhancement in the pro- 
tonic cosmic ray energy density is occurring due to local 
reacceleration. We estimate a luminosity of <~ 5.2 x 10 24 
erg s _1 for a cosmic ray enhancement of a factor 10, using 
equation (19) in Aharonian & Atoyan (1996). 



In the case of relativistic bremsstrahlung, which is 
produced in the electrostatic fields of local ions by the 
same electronic population that emits the observed syn- 
chrotron radiation, we have that the ratio of gamma-ray 
bremsstrahlung (F 7 ) to synchrotron radio (S syn ) fluxes is 
given by (see Pollock 1985): 



R 



F 7 (> 100 MeV) 



(5) 



~ 1.4 x 10~ 13 ( 
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The above expression has been derived for a canonical 
power-law electron spectrum of index p = 2 (N C (E) oc 
E~ p ). Using the observed radio flux at 5 GHz, and 
the same physical parameters than in the case of the 
7r°-decays, we obtain a bremsstrahlung contribution of 
~ 3.4 x 10 30 erg s _1 . This value is more than a thou- 
sand times smaller than the contribution from the inverse 
Compton emission, that clearly dominates the high-energy 
radiation from this region. 

We conclude that, according to the available radio ob- 
servations, a significant part of the observed gamma-ray 
flux from 3EG J2033+4118 might be produced in the col- 
liding winds region of Cyg OB2 No. 5. 

4-2. Terminal shock 

Particle acceleration and gamma-ray production in the 
shocked region at the boundary between a supersonic stel- 
lar wind and the surrounding medium has been discussed 
in detail by Casse & Paul (1980) and Volk & Forman 
(1982). White (1985) showed that energetic particles from 
the star can be injected into the acceleration region de- 
spite the strong adiabatic losses. At the terminal shock, 
the charged particles are accelerated to a power law dis- 
tribution with index 2. Observable gamma-ray radiation 
can be generated if there is a nearby source of photons 
to interact with the electrons or nearby clouds that might 
be "illuminated" by the relativistic protons. In the case 
of Cyg OB2 No. 5, if there is a contribution from the ter- 
minal shock to the total gamma-ray luminosity, it should 
come from the latter mechanism because there is no ev- 
idence of synchrotron emission from the terminal shock 
that reveals a particularly high energy density of rela- 
tivistic electrons. On the contrary, there is evidence of 
several compact clouds in the region. Dobashi et al. (1994, 
1996) have made an extensive study of molecular clouds 
in the Cygnus region. They have found 159 molecular 
clouds through large-scale 13 CO observations in the re- 
gion. The mass spectrum of these clouds is found to be 
well approximated by a power law dN/dM c \ oc M C T , for 
M c i > 100 M . In the neighborhood of Cyg OB2 No. 5, 
four clouds with estimated masses of > 230, > 6750, 1530, 
and < 240 Mq have been detected. The presence of many 
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protostellar candidates detected by IRAS also indicates a 
large number of small, dense clouds. It is not unreason- 
able, then, to assume that at least one small cloud of, 
say, 100 Mq could be illuminated by protons locally ac- 
celerated in the terminal shock of the wind, at a few pc 
from the primary star. If the local density of cosmic rays 
is enhanced in a factor 10 from the value observed in the 
Earth vicinity (e.g. Dermer 1986), we can expect an inte- 
grated gamma-ray luminosity of ~ 2.3 x 10 32 erg s _1 in 
the EGRET range, with the conservative assumption that 
only a small cloud of 100 Mq is overtaken. 

4-3. Base of the wind 

The winds of O and B stars are radiatively driven by ab- 
sorption in spectral lines and are prone to undergo insta- 
bilities that can grow up to form strong shocks at the base 
of the outflow (e.g. Lucy & White 1980, Lucy 1982, Lamers 
& Casinelli 1999). These shocks can accelerate particles 
to relativistic energies. If the particles are electrons, syn- 
crhotron and IC radiation can be produced. White (1985) 
have calculated the electron energy spectrum and the as- 
sociated synchrotron spectrum, and have shown that these 
spectra can be quite different from the naive pure power- 
law case. In particular, the wind synchrotron spectrum 
rises with frequency as S v oc iA 9 ^/ 6 , with (3 = 1.6 in the 
strong shock limit, for frequencies below a turnover v tl 
and is nearly flat (S v oc lni/ t ) from the turnover to the 
high energy cutoff. This theoretical synchrotron spectrum 
agrees notably well with the radio spectrum of the primary 
in Cyg OB2 No. 5, which presents a spectrum oc jA 54 be- 
tween 4.8 and 8.4 GHz (Waldron et al. 1998), that flattens 
at higher frequencies (Contreras et al. 1996). The value of 
ft, given by equation (45) of White (1985), is highly sen- 
sitive to small changes in the stellar parameters. For the 
values listed in Table 1 for the primary in Cyg OB2 No. 5 
we obtain t> t « 10 GHz, in good agreement with the radio 
data. 

The high-energy cutoff for charged particles acceler- 
ated in the vicinity of a star is (White 1985): 

which yields -E max ~ 2 GeV for the primary in Cyg OB2 
No. 5 (i.e. Lorentz factors 7 ~ 4 x 10 3 ). This upper limit, 
imposed by geometric factors, is valid for both electrons 
and ions. However, in the case of electrons, radiative losses 
mainly due to IC scattering impose additional constraints 
(Chen & White 1991), and hence the effective high-energy 
cutoff should occur at 7 ^ 10 3 for them. Consequently, 
electron-induced gamma-rays in the base of the wind have 
energies Ss 100 MeV and are below EGRET energy thresh- 
old, as noticed by Chen & White (1991). 

Relativistic ions, on the contrary, can produce sub- 
stantial gamma-ray luminosities through hadron interac- 
tions with the thermal wind nuclei and the subsequent n°- 



decays (White & Chen 1992). The resulting high-energy 
emission has a peak at <~ 70 MeV and then can sig- 
nificantly contribute to the gamma-ray flux detected by 
EGRET. 

The calculation of the ir° gamma-ray luminosity of the 
primary star in Cyg OB2 No. 5 has already been carried 
out by White & Chen (1992), using the empirical formula 
for the 7T° production spectrum given by Stephens & Bad- 
hwar (1981). The expected luminosity is ~ 5 x 10 34 erg 
s _1 , similar to the IC luminosity of the colliding winds 
region. 

4-4- Summary of results 

In Table 2 we present a summary of the different contribu- 
tions to the total gamma-ray luminosity of Cyg OB2 No. 
5 in the EGRET energy range. The main contributions 
are IC emission from the colliding winds region and ir°- 
decay emission from the inner region of the primary stellar 
wind. These mechanisms can provide about a half of the 
inferred luminosity of the source 3EG J2033+4118 (i.e. 
<~ 1.3 x 10 35 erg s _1 ). The observed gamma-ray flux could 
be explained if either (1) the mass illuminated by relativis- 
tic protons accelerated at the terminal shock of the wind 
is a factor ~ 50 larger than what we have assumed (i.e. if 
it is ^5000 M Q ), (2) there are other gamma-ray sources 
within the EGRET 95 % probability location contour or 
(3) there is an important contribution from the colliding 
winds region within the primary (i.e. in the close binary). 
This latter possibility cannot be quantitatively discussed 
because of the absence of observational constraints on the 
efficiency of the particle acceleration at the unresolved in- 
teracting region. The second possiblity, namely that the 
EGRET detection could actually be a composed source, 
will be briefly discussed below. 

4-5. Other potential gamma-ray sources in the field 

There are, in addition to Cyg OB2 No. 5, other massive 
stars of the same association within the positional error 
box of 3EG J2033+4118, like stars No. 7, 8A, 9, and 11. 
These stars have similar parameters to the primary in Cyg 
OB2 No. 5 and are potential sources of 7r°-decay gamma- 
rays. In particular, Cyg OB2 N0.8A and No. 9 are non- 
thermal radio emitters (Contreras et al. 1996, Waldron et 
al. 1998). White & Chen (1992) have estimated the indi- 
vidual gamma-ray luminosities at E > 100 MeV of these 
stars in ~ 5 x 10 34 erg s _1 . The combined effect of these 
stars, plus the emission from the colliding winds region 
in Cyg OB2 No. 5, can easily explain the flux detected 
by EGRET. The absence of any other strong radio source 
in the field also supports this conjecture, which was ad- 
vanced by White & Cheng (1992) prior to the detection 
of non-thermal radio emission from Cyg OB2 No. 5 (see 
also Chen et al. 1996). 
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5. Further comments 

5.1. High-energy tail of the gamma-ray emision 

In the colliding winds region of Cyg OB2 No. 5 particles 
are accelerated at the shock front up to a high-energy 
cutoff that can be estimated assuming that the diffusion 
length is equal to the gyral radius. This yields, 

£ max = 1.3 x 10 3 GeV r-j-\ r 2 , (7) 

where v s is the shock velocity in km s -1 . If v s = Voo,i and 
B = B s , we obtain £ max ~ 0.7 TcV. The IC high-energy 
emission, then, extends up to energies ~ 20 TeV. The 
expected gamma-ray luminosity at E > 1 TcV is 5 x 10 30 
erg s _1 , which, at a distance of 1.8 kpc gives a photon 
flux of 1.4 x 10~ n cm _2 s _1 . This flux could be detected by 
planned gamma-ray ground arrays like VERITAS (Weekes 
et al. 2000). 

5.2. A few comments on the soft gamma-ray emission 

The colliding winds region can be an important source of 
soft gamma-ray photons. The IC flux density around 10 
MeV is 8 x 10~ 5 ph cm~ 2 s _1 and could be detected, in 
principle, by the Imager on Board INTEGRAL Satellite 
(IBIS), to be launched by the European Space Agency in 
April 2002 (Winkler & Hcrmscn 2000). IBIS is a coded 
mask telescope designed to produce sky images in the 20 
KcV - 10 MeV band with an angular resolution of 12 
arcminutes. It will also have a good spectral capability. 
The point source location accuracy for a relatively strong 
source will be around 1 arcminute. Consequently, IBIS will 
provide excellent diagnosis to test the models for gamma- 
ray production in Cyg OB2 No. 5 discussed in this paper. 
In particular, it will allow to measure the soft gamma-ray 
spectrum of the colliding winds and differentiate the emis- 
sion of Cyg OB2 No. 5 from that produced in the neigh- 
boring stars, which are located at an angular distance > 9 
arcminutes. The main problem at these energies is the 
subtraction of the background galactic emission; notwith- 
standing, sufficiently accurate models of the diffuse radi- 
ation fields at low latitudes will be soon available (e.g. 
Blocmen et al. 2000). IC emission from the single stars 
should also be important at MeV energies. The hadronic 
gamma-ray component, however, should be observed with 
future higher energy telescopes, like the planned GLAST 
satellite, to be launched by NASA towards 2005. 

6. Conclusions 

We have presented an analysis of the possible produc- 
tion of gamma-ray emission in the stellar system Cyg 
OB2 No. 5 taken into account the constraints imposed 
by recent high-resolution radio observations. We have es- 
timated the gamma-ray luminosity generated by different 



types of mechanisms and compared them with the lumi- 
nosity inferred from EGRET observations. We have found 
that about a half of the flux of the gamma-ray source 
3EG J2033+4118 could arise from the colliding winds re- 
gion (through IC interactions) and the base of the wind 
(through proton- proton interactions). It is possible, also, 
that the remaining flux comes from other nearby early- 
type stars of the same OB association within the posi- 
tional error box of the EGRET detection. We expect that 
forthcoming gamma-ray missions with good imaging capa- 
bilities, like the imminent ESA's INTEGRAL spacecraft, 
can clarify this point. 
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Table 1. Stellar parameters 



Star Sp. Class. logT eff R* log(L/L ) B t M V ro t 
(R©) (G) (lQ- 5 M /yr) (km/s) (km/s) 

Primary 07 Ianfp (a) 4.6 (c) 34 (d) 5.984 (c) 2 2.5 (f) 2200 (d) 180 (d) 

Secondary BO V (b) 4.5 (c) 8.3 (c) 4.881 (c) 1 0.2 (f) 1000 (b) 150 (o) 

(a): Walborn 1973; (b): Contreras et al. 1997; (c) Vacca et al. 1996; (d): Bieging et al. 1989; 
(e): Uesugi & Fukuda 1982; (f): See text. 



Table 2. Gamma-ray production in Cygnus OB2 No. 5 for EGRET energy range 
(100 McV < E < 20 GcV) 



Region 


Mechanism 


Expected luminosity 
(erg/s) 


Observed luminosity 
(erg/s) 




Inverse Compton scattering 


~8x 10 34 




Winds collision 


Relativistic bremsstrahlung 


~ 3.4 x 10 30 






Neutral pion decay 


~ 5.2 x 10 24 


~ 2.4 x 10 35 


Terminal shock 


Neutral pion decay 


~ 2.3 x 10 32 




Base of the wind 


Inverse Compton scattering 








Neutral pion decay 1 


~ 5 x 10 34 





1: White & Chen 1992 



